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The spectrum of supersymmetric theories with R-parity violation are much more weakly con- 
strained than that of supersymmetric theories with a stable neutralino. We investigate the signatures 
of supersymmetry at the LHCb experiment in the region of parameter space where the neutralino 
decay leaves a displaced vertex. We find sensitivity to squark production up to squark masses of 
order 1 TeV. We note that if the Higgs decays to neutralinos in this scenario, LHCb should see the 
lightest Higgs boson before ATLAS and CMS. 



I. INTRODUCTION 

Supersymmetry, and more specifically the minimal su- 
persymmetric standard model (MSSM) [l| is a possible 
solution to the gauge-hierarchy problem and a favorite 
model governing physics above 100 GeV. Quantum cor- 
rections to the electroweak breaking scale are propor- 
tional to the superpartner masses, and thus one expects 
the MSSM spectrum to lie around the Z mass. Any 
significant deviation thereof necessitates a fine-tuning of 
parameters. 

Minimal Supergravity (mSUGRA) is the most studied 
realization of the MSSM. It assigns universal masses to 
all scalars and to all gauginos at the scale 10 16 GeV. Ex- 
perimental bounds on the mSUGRA spectrum demand 
that the model is tuned to the per cent level [|J • Squarks 
and gluinos have lower bounds around 300-400 GeV 
well above the Z mass. The corrections necessary to gen- 
erate a Higgs mass consistent with direct LEP bounds Q 
requires even more fine-tuning. The tight constraints on 
mSUGRA compel us to study more natural - i.e. less 
fine-tuned models of supersymmetry. 

The most general superpotential with the MSSM field 
content includes lepton and baryon number violating 
terms [f| 

XiikLiLjEl + \' ijk LiQjD%. + \': jk UfDpi, (1) 

where ijk are flavor indices. Bounds on proton decay 
severely constrain the combination of baryon and lepton 
number violation. Separately, however, they are much 
more weakly constrained. Throughout the following we 
restrict our attention to the baryon number violating 
operators, the A terms. This choice is motivated by 
the interesting and challenging phenomenology it pro- 
duces. The bounds for squark masses in this scenario 
are typically below 100 GeV @; some particles, such as 
the gluino and lightest sbottom, do not have published 
bounds above ~ 10 GeV in regions of parameter space 
Q. Moreover, the bound on the Higgs can be below the 
Z mass when decays to neutralinos are kinematically al- 
lowed @|. 

The phenomenologically interesting feature of this 
model is that the lightest superpartner - taken to be a 
neutralino - is unstable. It decays to three quarks. With 
regard to supersymmetry at the LHC, the signals are 
changed significantly - missing energy signals are largely 



absent, and the number of isolated leptons is reduced 
due to increased soft jet production flfj| |. The ATLAS 
and CMS experiments 0] will have weakened sensitivity 
to this scenario because their triggers are designed to ex- 
ploit missing transverse energy and isolated leptons. In 
the case of squark or gluino production the associated 
jets should pass the triggers at ATLAS and CMS, how- 
ever hard jets are typically pre-scaled by a large factor 
that would significantly reduce the effective luminosity 
[111 ]. Even if the trigger issue is solved, it is not clear 
that pure multi-jet events coming from this new physics 
can be seen above the (unknown) QCD background. Yet 
more worrisome is the Higgs decay. For decays to a final 
state of six soft jets, no obvious search strategy presents 
itself, while the standard searches are made more difficult 
with the reduced branching ratios. 

The lightest neutralino has a macroscopic decay length 
in broad regions of parameter space. While neither AT- 
LAS nor CMS currently employ a displaced vertex trig- 
ger, LHCb 12] is designed to trigger on and reconstruct 
such events. LHCb operates at a center of mass energy 
equal to that of ATLAS and CMS (14 TeV). However, 
its luminosity is limited to 2 fb _1 per year and it cov- 
ers only the forward region. The experiment is designed 
to make measurements of rare 6-hadron decays by relying 
on their ability to precisely reconstruct displaced vertices. 
The lower luminosity limits the average number of inter- 
actions per bunch crossing to < 1, which allows for more 
precise vertexing. This makes it an ideal experiment to 
search for our signal. 

The purpose of this article is to show quantitatively 
that the LHCb experiment should have significant reach 
in the parameter space of this class of supersymmetric 
models. In parts of parameter space, it may be able to 
show the first direct evidence of the lightest neutralino 
and the lightest Higgs boson. Below we present our esti- 
mates of the LHCb's physics reach with regard to squark 
and Higgs production. Our work shows that the signal 
events easily pass the lowest level LHCb triggers. We 
suggest a modified high level trigger to increase the ef- 
ficiency with which the signal is written to tape. While 
computational limits prohibit our complete understand- 
ing of the leading order QCD background, we argue that 
for some parts of parameter space the signal will domi- 
nate the background. 
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II. NEUTRALINO DECAY VIA BARYON 
NUMBER VIOLATION 

The baryon number violating operators in Eq. JT]) in- 
volve nine complex couplings (because j and k are anti- 
symmetric). When the neutralino decays, it does so via 
the coupling X". k into up- type quark i and down- typed 
quarks j and k through an off-shell squark. A reason- 
able, theoretically motivated parameterization for these 
couplings based on a spurion analysis of flavor breaking 
in the standard model is 11311 



A ijk 
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mimjirik 
v 3 sin /3 cos 2 /3 ' 
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where the m,, etc., are quark masses jl4| . and v sin /3 and 
v cos /3 are the vacuum expectation values of the up- type 
and down- type Higgses respectively, with v — 17 A GeV. 
We shall use this parameterization and take tan/3 = 1 
since any difference can be absorbed into An. Note, the 
A323 coupling dominates. The dominant decay mode of 
the neutralino will be \ ~~ * ^bs : unless the neutralino is 
lighter than the top in which case x ~* c ^ s dominates. 
In either case, neutralino decays are dominated by heavy 
flavors, and should contain additional displaced vertices. 
We will not utilize this additional handle on the signal, 
though it may prove to be a useful part of the full exper- 
imental analysis. 

The strongest current bounds on the magnitude of A" 
couplings are from baryon number violating processes, 
namely neutron-antineutron oscillations and double nu- 
cleon decay in, for example, oxygen nuclei [HI]. Such 
bounds allow A ' ~ 0(1) within QCD uncertainties. If 
the A" have arbitrary complex phases, they can con- 
tribute to direct CP violation in Kaon decays and to 
K — K mixing. The strongest bound in this case is the 
limit 2^3.^323*) < 10 -8 [lq , which implies a bound on 
our universal parameter Aq < 1/20 if all phase differences 
are order unity and squark masses are 100 GeV. There 
are no significant bounds on the individual A 223 and A 3 ' 2 3 
couplings. For a complete review, see [l7| . 

The proper lifetime of the neutralino depends on the R- 
parity violating couplings, the neutralino mass ro x , and 
the squark masses m,q. With the simplifying assumptions 
of a universal squark mass at low energies and a single 
dominant R-parity violating coupling (as in our parame- 
terization), the proper lifetime is 
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where |?7 2 i| is an element of the neutralino rotation ma- 
trix (see Q). We have neglected Yukawa couplings, 
QCD corrections and phase-space corrections (taking fi- 
nal state particles as massless). These are good approx- 
imations in the two cases we study. Yukawa couplings 
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FIG. 1: A qualitative picture of a neutralino decay off the 
beamline. The decay is of a neutral particle (no track) into a 
large multiplicity of tracks. The invariant mass of the tracks 
should be significantly larger than those from a &-hadron de- 
cay. The quantities r min and z min are defined in the text. 



are relevant to the extent that the lightest neutralino is 
partially higgsino. For Higgs production and decay, the 
neutralino is much lighter than the top, and for the de- 
cay to neutralinos to dominate, it turns out tan/3 should 
not be too large Q, and therefore all relevant Yukawa 
couplings are small. In the case of squark production, 
we will look only at the 'pure bino' limit (making the 
higgsinos and winos heavy), and thus we can ignore the 
Yukawas entirely. In the former case, |J7 2 i| is less than 
and of order unity. In the latter case (pure bino limit), 

1^21 1 = 1- 



III. SIGNALS AND BACKGROUNDS AT LHCB 

Here we list the expected signals and backgrounds, pro- 
posed triggers and signal efficiencies, and offline discrim- 
inants: 

• In Figure[T]we show pictorially a macroscopic decay 
of a neutralino off the beam line. The production 
signals we study are: 

— Squark production with q — > qxo and \a 
qqq. 

— Higgs production with h — > xaXo an d Xo ~^ 
qqq. 

They are generated, including showering, with 
Pythia v6.400 [18]. For squark production we use 
the following parameters: 

— The ratio of couplings (A 223 /A 323 ) = (1/20). 

— A scan of m s from 100-1000 GeV in 100 GeV 
steps. 

— A scan of three different bino masses: Mi = 
40, 100, 200 GeV and three different coupling 
values A 223 = 10~ 3 , 10~ 4 , and 10" 5 . 

— Mi = M 3 = (i = 1.2 TeV while all other pa- 
rameters are set to the Pythia default values. 
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For Higgs production we do the same scan of Mi 
parameters and use tan/3 = 5, M 2 — 250 GeV, 
fx = 120 GeV, TTiq = A t = 1 TeV and A 2 ' 23 = 10~ 2 , 
with other soft terms at default values and other 
A" couplings set to zero. 

• The background is taken to be multiple b produc- 
tion. We use Pythia to simulate bb events. Mad- 
graph v4.1.19 [19| is used to compute matrix ele- 
ments of gb — > bbb, gg — ► bbbb, and gg —> bbcc while 
Pythia is used to shower these events. 

• We find that the following cuts and triggers dis- 
criminate the signal and background: 

— Requiring a displaced vertex with 300 /im < 
z < 0.4 m and r > 60 /im, where z and r 
are the horizontal and perpendicular distance 
from the interaction point. 

— Requiring at least 5 tracks from the displaced 
vertex. 

— Requiring at least two tracks with pt > 1 GeV 
and a two-dimensional impact parameter of 
0.07 mm < bjp < 15.0 mm 

• For offline discrimination we use invariant mass dis- 
tributions of displaced vertices. 

We now describe the relevant aspects of the LHCb exper- 
iment [l2j and explain in detail the motivation for and 
expected results of these cuts. LHCb is asymmetric in 
theta (polar angle) acceptance. The horizontal accep- 
tance is 15 mrad < 8 < 300 mrad while the maximum 
vertical acceptance is 250 mrad. For simplicity we restrict 
our analysis to the region 15 mrad < 8 < 250 mrad. Of- 
fline reconstruction of the primary vertex is expected to 
have a resolution of < 50 /mi along the beam line and < 
10 /im perpendicular to it . The typical z resolution of a 
secondary vertex is ~ 200 /im. Transverse resolution is 
Pt dependent, and is ~ 20 /im for 1 GeV pr track. We 
assume that vertexing may be done up to 0.4 m along the 
beamline which corresponds to half of the Vertex Locator 
length 20], We set the resolution of a displaced vertex 
to be a cylinder of 200 /im in z and 30 /im in r. This 
means that if a second vertex lies outside this cylinder 
then it can be distinguished, otherwise it cannot. We de- 
note these lengths as 8z and Sr. The required minimum 
distances from the primary vertex as described above in z 
and r are denoted z m i n = 300 /im and r rn i„ = 60 /mi. No 
detector effects beyond vertex resolution are considered. 



A. Squark Production 

All superpartners produced at the LHC cascade to the 
lightest neutralino (direct decays via R-parity violation 
are suppressed by the small A" coupling). For simplicity 
we look at squark pair production where each squark 
decays to a quark and the lightest neutralino. The goal 
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FIG. 2: Efficiency of xo acceptance with respect to impact 
parameter requirements. We require 5 tracks, 2 of which with 
Pt > 1 GeV. The lower points show 0.15 mm < bip < 3.0 
mm and the upper points show 0.07 mm < bip < 15.0 mm. 
The parameters are set at A 223 = 10" 4 , Mi = 200 GeV, and 
mq = 400 GeV. 



of this search is to see one of these neutralinos in the 
LHCb acceptance. The signal we look for is a displaced 
vertex with a larger track multiplicity than a typical b 
decay and an invariant mass of the tracks larger than the 
b mass. The decay length may or may not be similar to a 
typical b hadron, so we do not use this as a distinguishing 
feature. 

LHCb's Level (L0) trigger is designed to reject mul- 
tiple primary vertices ('pile up' events) and events with 
large numbers of tracks ('busy' events). L0 reduces the 
data rate from ~ 12 MHz to - 1 MHz. We find that 95% 
of our squark production signal (one neutralino leaving at 
least 5 tracks in the detector acceptance) passes L0. The 
High Level Triggers (HLT) are responsible for reducing 
the rate to 2 kHz, the read out rate. One component of 
the HLT is 2D track reconstruction searching for tracks 
with high pt (> 1 GeV) and large impact parameter, 
0.15 mm < bjp < 3.0 mm, tracks. Our signal generically 
produces more high pt tracks than the background be- 
cause of the neutralino's greater mass and because it is 
the product of a heavy particle decay. However, we find 
that the proposed impact parameter window results in 
signal efficiencies below 10% for decay lengths inconsis- 
tent with that of a b. The signal efficiency is increased to 
above 50% in most parts of parameter space if the impact 
parameter window is widened to 0.07 mm < bjp < 15.0 
mm. Figure shows the efficiencies of the two ranges 
for a particular point in parameter space. 

Extending the impact parameter range to a lower value 
of 50 /xm is suggested in the context of LHCb upgrades 
[2l| . The feasibility of extending the range to large val- 
ues in unknown and requires a detector simulation. The 
naive background for large impact parameter tracks is 




FIG. 3: Number of expected \o events from squark produc- 
tion vs. squark mass. At least 5 tracks with 2 having more 
than 1 GeV of pr and 0.07 mm < bip < 15.0 mm are required. 
The coupling is A223 = 10~ 4 and the neutralino masses are 
computed by Pythia using the parameters set at the begin- 
ning of the section. 
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FIG. 4: The same as Fig [3] but with A 223 = 10" 



strange decays. This is because T strange c ~ centimeters. 
If this is the case, it seems plausible that requiring mul- 
tiple tracks with high px can significantly reduce this 
background. Henceforth we optimistically assume that 
the range can be extended. 

The expected event rate of neutralinos from squark 
production that pass our cuts are shown in Figures [3] 
and |U 

Figure O shows the signal invariant mass distributions 
of tracks for events that pass our cuts. As expected, these 
distributions reach far beyond the mass of a 6-hadron, 
and thus we need not worry about isolated 6 decays as 
a background. We find rather that the most important 
background (offline) consists of 6-hadron events in which 



FIG. 5: xo Invariant Mass from squark production. All points 
are A223 = 10~~ J with the same requirements as Figure |3] 
Red: Mi = 40 GeV, m§ = 100 GeV. Blue: Mi = 100 GeV, 
m§ = 200 GeV. Green: Mi = 200 GeV, ma = 400 GeV. 



another particle, produced at the primary vertex, decays 
near enough to the 6 that the vertices cannot be resolved 
individually by the detector. We refer to these as over- 
lapping events. 

We now give a rough quantitative estimate of the back- 
ground. Our region of interest is for track invariant 
masses above 2ms ~ 12 GeV. We see from Figure [5] 
this region has a significant overlap with our signal. Our 
limited computing power only allows us to simulate 10 -5 
years of background, in which we find no events which 
pass our track cuts and have an invariant mass above 
12 GeV (see Figure IB)). To better understand the back- 
ground, we also look at the invariant mass of all decay 
products (charged and uncharged) from these overlap- 
ping events. Using this information, we are able to define 
cuts which should in principle reduce our background to 
less than 1000 events per year. Below, we describe how 
we come to this estimate. 

The expression for the invariant mass of two particles 

is 



M 2 = m\ 



2(ExE 2 -p Zl Pz 2 -pTrP^cosAc/)). (4) 



where mj, E%, p Zi , and pn are the mass, energy, z- 
momentum and transverse momentum respectively of the 
ith particle, and A<j> — 4>\ — 4>2 is the difference in the 
azimuthal angle of the two particle momentum vectors. 
There are two overlap cases: 26s and a b plus a non-6. We 
discuss the 26 case. It is clear that the non-relativistic 
limit cannot produce M 2 3> (2m&) 2 . The relativistic 
limit reduces (J4j) to 
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where Cj = cos 6*;, the cosine of the polar angle of the 
momentum vector of the ith particle, and similarly, Sj = 
sin^j. (Note, the region of parameter space where one 
6 is non-relativistic is a special case of what we discuss 
below) . 

Examining the cross term we see that there are two 
interesting cases: p Zl ~ p Z2 and p Zi 3> p z ■ The for- 
mer case requires a large difference in polar or azimuthal 
angles to generate a large cross term. Maximizing this 
difference (for example, in the polar angle) while demand- 
ing a large cross term and using px = p z tan 9 leads to a 
minimum px for the 6— hadrons. Furthermore a large 9 
difference with a small transverse distance (making the 
vertices unresolved) requires the vertices to be as near 
the primary vertex as possible. These considerations sig- 
nificantly suppress the number of overlapping events. As 
an illustration we take 8r = 30/xm, the closer b a trans- 
verse distance of 30 /im from the z axis, and the vertices 
a distance of 360 /im from the primary vertex along the z 
direction. Demanding the cross term give f 6 rn^ (to get 
an invariant mass of all decay products, not just tracks, 
of just over 20 GeV) we find that p z > 160 GeV. This 
corresponds to pr > 20 GeV for the softer b, a require- 
ment which suppresses the cross section by better than 
10 -5 and makes this parameter range irrelevant. 

Conversely, the case in which p Zi 3> p Zj is important 
even when the 9 difference is small. The non-relativistic 
corrections - the last term in Eq. ([5]) - dominate when 
the difference in angles vanish. Generating a cross term 
of 16m^ requires a ratio of 16:1 between the p z s. The 
softer b (call it 6i) decays dominantly at a length L < 
(p^/m^TbC. Now the harder b has pr 2 — p Z2 9 ~ \Qp Zl 9. 
Using the requirement that 9 > ( so th e displaced 

vertex satisfies our r m i n cut), and plugging the values 
of r m in and Tf,c leads to the requirement pr 2 > 2mb. 
To estimate our background, we create a sample of 10 -5 
years of 26 production using Pythia (roughly 10 7 events 
at leading order) requiring both 6s to decay within the 
acceptance of our detector and to pass our r m in and z m in 
requirements. We then count the number of events that 
satisfy (p Z2 /p Zl ) > 16 and pr 2 > 10 GeV. The fraction of 
our sample which passes these cuts is one part in 2 x 10 3 . 
Then we take the same sample without the momentum 
requirements and find the number of overlap events to be 
69 - or scaled up, roughly 7 x 10 6 per year. If we take the 
distribution of momenta among these events to be flat 
(overly conservative), we can simply take a product of 
the two suppressions and estimate the number of events 
which have the potential to have a large enough invariant 
mass. Our estimate is N < 7 x 10 6 x 5 x 10~ 4 = 3, 500. If 
we include the fact that 6's with very different momenta 
will have very different decay lengths, we find another 
suppression of a factor of nearly an order of magnitude 
and thus expect a background to be at most on the order 
of hundreds of events. 

In addition to the 66 background, there are overlap 
events generated in, for example, the 36 background. We 
find no events where two 6s overlap giving a large invari- 
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FIG. 6: Background invariant mass of overlapping events in 
2b and 36 production. Objects are defined as all decay prod- 
ucts that deposit energy in the detector (i.e. everything but 
v's). We demand more than 5 objects or more than 5 tracks 
from a vertex. No cuts are made on impact parameter nor 
track/object pr- All large mass events in red are due to b 
decay overlapping with a non— b decay. 

ant mass, and using similar arguments to those above 
find that a full year should produce at most as many 
events as in the 26 sample. However, we do find large 
invariant mass events in this sample which involve the 
overlap of a 6 and a strange hadron decay. These events 
do not pass either the 5 track cut or the r m i n require- 
ment. Imposing both should in principle severely limit 
or eliminate events of this type, but unfortunately it is 
difficult to estimate. We will simply assume they can be 
removed by these or similar cuts. In Figure[S]we plot the 
invariant mass of overlapping vertices in the 36 sample. 
We include in the plot the invariant mass of all decay 
products to see the large invariant mass events. All of 
the events with invariant masses larger than 10 GeV are 
due to 6-non-6 overlapping events. 

The other simulated backgrounds^ produce an overlap- 
ping event fraction about that of 66, and they are cross 
section suppressed by more than an order of magnitude. 

We now estimate the parameter space that can be ex- 
plored by LHCb. We assume a background of 400 events 
above 12 GeV of track invariant mass. Significance at 
the level of -7= > 5 requires > 100 signal events above 
12 GeV. The regions of parameter space which exceed 
this event rate after cuts are: 

• a;' 23 - 10- 3 : 

- M X0 = 38: 200 GeV < m 5 < 600 GeV. 

• A; 2 3 - f0- 4 : 

- M X0 = 38: 100 GeV < m 9 ~ < 400 GeV. 

- M X0 = 98: 200 GeV < m s < 700 GeV. 
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FIG. 7: Number of Expected xo in the acceptance from Higgs 
production and decay at 2 fb _1 . We require 5 tracks, 2 of 
which with more than 1 GeV of pr and 0.07 < bip < 15.0 
mm. This point is tan/3 = 5, M 2 = 250 GeV, fi = 120 GeV, 
rriq = A t = 1 TeV, and all other parameters at Pythia default 
values. The mass of the Higgs at this point is van ~ 115 GeV. 

• a;' 23 = nr 5 : 

- M xo = 38: 100 GeV <m^< 200 GeV. 

- M Xo = 98: 200 GeV <m^< 400 GeV. 

- M Xa = 198: 300 GeV < m 9 < 700 GeV. 

B. Higgs Production 

There exists an interesting region of parameter space 
in which the Higgs dominantly decays to neutralinos. 
The signal invariant mass distribution becomes more 
background-like in this parameter space as m Xo becomes 
lighter. In fact the distribution is almost indistinguish- 
able from the overlap background in Figure [5] for m Xo < 
20 GeV, in part because the neutralino vertex loses some 
of its track invariant mass through its decay to a & whose 
decay products often reconstruct at a different point. 

However, a distinguishing characteristic of the signal 
is the possibility of both neutralinos being in the accep- 
tance. An exemplary point is shown in Figure [7] where 
we use the leading order inclusive Higgs production cross 
section of Oh = 20 pb [22]. We see that at this point 
if m Xo > 25 GeV then there are a significant number of 
Higgs decays that deposit both neutralinos in the detec- 
tor. The expected background for this signal is two sets 
of overlapping decays which should be negligiablc. 

The rapid fall off of the acceptance distribution at 
small m Xo is primarily due to long decay length at these 



particular points. The decay length of a 20 GeV neu- 
tralino for the prescribed parameters is L ~ 0.2 m. Given 
that such a neutralino is boosted by at least a factor 3 
due to the Higgs decay, we would not expect both neu- 
tralinos to be in the acceptance (defined with an upper z 
limit of 0.4 m) whereas it is not surprising that one of the 
decays is in the detector. Reducing m Xo by a factor of 
2 leads to 32 fold increase, plus a significant phase-space 
increase of L - hence the sparseness of accepted events 
at small rn Vn . 



IV. DISCUSSION 

The addition of baryon number violating operators 
to the MSSM superpotential allows for a more natu- 
ral model of supersymmetry while also producing phe- 
nomenology that may pose difficulties for ATLAS and 
CMS. The central phenomenological signature is dis- 
placed vertices for which the LHCb is well suited to ob- 
serve and reconstruct. Displaced vertices of b— hadron 
decays is a potentially enormous background. Despite 
computational limitation that forbid a full simulation of 
the background we can argue that it is plausible for large 
portions of parameter space to be explored. This is a 
consequence of the large invariant mass distribution of 
our signal, that necessitates coincident background de- 
cays. Thus we estimate that LHCb could rule out a sig- 
nificant portion of parameter space. However, only the 
most naive detector issues have been considered and a 
full detector simulation is needed to understand the de- 
tector's true reach. 

In addition to R-parity violating supersymmetry, other 
versions of supersymmetry may also contain displaced 
vertices. This includes parts of parameter space with 
near degeneracies between the LSP and NLSP which 
could occur between, for example, a stau and a neu- 
tralino, or between neutralinos in theories with an added 
singlet field. Finally, so-called 'hidden valley' models [23[ 
also give rise to non-standard displaced vertices and have 
been suggested as good candidates for LHCb physics. A 
dedicated search at LHCb may provide the first discovery 
of new physics at the LHC. 
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